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ABSTRACT 

We describe a mechanical equation of state for NH3, based on shock wave 
measurements for hquid ammonia. The shock measurements, for an initial tem- 
perature of 203 K, extended to 1.54 g/cm^ and 38.6 GPa. The shock and particle 
speeds were fitted well with a straight line, so extrapolations to higher com- 
pressions are numerically stable, but the accuracy is undetermined outside the 
range of the data. The isentrope through the same initial state was estimated, 
along with its sensitivity to the Griineisen parameter. Mass-radius relations 
were calculated for self-gravitating bodies of pure ammonia, and for differen- 
tiated ammonia-rock bodies. The relations were insensitive to variations in the 
Griineisen parameter, indicating that they should be accurate for studies of plan- 
etary structure. 

Subject headings: ammonia, shock, equation of state, planetary structure 
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Introduction 



Ammonia, NH3, is a com mon molecule in ice giant planets, which appear to be found 



widely throughout the galaxy ( ISchneider 



201l[ ). The equation of state (EOS) of ammonia 



is therefore important for our understanding of planetary structures and their evolution, 
potentially to pressures of order 1 TPa for the base of the ice-rock interface in icy exoplanets. 
An accurate EOS for ammonia is also needed for studies of hypervelocity impacts, such 
as meteoroid collisions with ice giants. Furthermore, ammonia is a simple prototype for 
bonds occurring in chemical explosives, for which densities from up to around twice that of 
zero-pressure solids are of interest for shock initiation and detonation. 

Although shock compressi on experimen ts have been performed on ammonia to pressures 



of several tens of gigapascals ([Marsh 



is SESAME 5520 



( iHaar fc Gallagher 



lolian 



1980l ). the only equation of state readily available 



1984 ). based on National Bureau of Standards gas phase data 



19781 ). and is tabulated to a maximum density of 0.765g/cm'^, which is 
barely greater than the zero-pressure density for liquid ammonia. Quasistatic compression 
experiments hav e been performed in which the density and sound speed were measured 



along isotherms (lAbramson 



2008 



Lieta] 



20091), but the highest pressures reported have 



reached only a few gigapascals. 



2. Empirical Griineisen equation of state 

Shock experiments have been reported previously in which the shock and particle 
speeds Us and Up were measu red for a ran ge of shock pressures, for liquid ammonia at an 



initial temperature of 203 K (iMarsh 



19801 ) (Fig. [T]). The uncertainties in Ug and Up were 



approximately 1%. These data can be fitted by a straight line fit 



Co + SiUp 



(1) 
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with 



Co = 2.00 ±0.13 (6.7%) 
si = 1.511 ±0.039 (2.6%), 

where the standard errors shown are from the residual fitting error, neglecting the 
uncertainty in measurement. The experimental measurements exhibited a slightly curved 
trend, but the number of points was not great enoug h to justify a higher-order fit. Solving 



the Rankine-Hugoniot equations for a steady shock (jZel'dovich fc Raizer 



19661 ) using the 



fitted parameters rather than the individual shock measurements, the highest shock pressure 
was 38.6 GPa, giving a mass densi ty of 1.54 g/cm^. The obs erved sound speed in liquid 



ammonia at 203 K is 1.9535 km/s (IBowen fc Thompson! 119681 ). which is consistent with the 



extrapolated Hugoniot data. (Fig. [T]) 

Various universal EOS have been proposed for differ ent classes of materia l. It is 



Woolfolk et all fll973[ ). whose 



interesting to compare with the 'universal liquid EOS' of 
only material-specific parameter is the sound speed at zero pressure. This EOS does not 
reproduce the shock data for ammonia, which is softer and more linear than the universal 
EOS (Fig.HD. 

The fit to the shock Hugoniot can b e used to predict the mechanical equation of state, 



using the Hugoniot as a reference curve (jMcQueen et al 



1970) 



p{p,e) 

Pr{p) 

er{p) 



Pr{p) ± r(p) [e - er{p)] 

Clprpjp - Pr) 

\p^sx{p- p,)f 

eo ± ^VAP) 

2 \pr p 



(2) 
(3) 

(4) 



where p^ is the initial density on the reference curve, and p.r(p) was derived for zero initial 
pressure, as here. Other experiments are required to determine r(yo), such as sound speed 
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Fig. 1. — Principal shock Hugoniot of liquid ammonia (initial temperature 203 K): experi- 
mental measurements and least-squares fit. The point at zero particle speed is the observed 
sound speed, which was not included in the fit. The curve labelled WCS is the Woolfolk- 
Cowperthwaite-Shaw universal liquid equation of state, whose sole fitting parameter is the 
sound speed at zero pressure. 
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measurements on the Hugoniot, a shock Hugoniot from a different initial state, or ramp 
compression. However, F can be est imated from the slope of the shock Hugoniot as 2si — 1, 



which is accurate for cubic crystals (ISkidmore 



19651). Thus pr = 0.725 g/cm^ and T ~ 2.022. 



Given the mechanical EOS, the ise ntrope thro ugh any state can be calculated by 



integrating the —p dv work numerically (jSwift 



20081 ). Isentropes calculated from Griinseisen 



EOS fitted to shock data typically behave unphysically at high compression, where the 
assumption that the Griineisen parameter is a function of density only breaks down. For the 
ammonia fit, the breakdown occurred at a mass density of 2.145g/cm'^. The isentrope was 
well-behaved to several terapascals, though its accuracy was undetermined. The isentrope 
should be reasonably accurate at least up to the peak compression in the shock data, 
which equates to 22.1 GPa on the isentrope. To investigate the sensitivity to the Griineisen 
parameter, isentropes were calculated for the nominal value above, and for values 10% lower 
and higher. With this variation in F, the pressure varied by 10% at 20 GPa, rising to 25% 
at 500 GPa. (Fig. |2]) 



3. Mass-radius relationships 



Mass-radius relationships were calculated using the deduced EOS, for a self-gravitating 



body comprising pure ammonia and also for differentiated bodies consisting o: 



and an ammonia mantle, using the numerical methods described previously (ISwift et al 



a rocky co re 



201ll ). Separate mass-radius curves were constructed for the nominal and perturbed values 
of the Griineisen parameter. In all cases, the temperature at the surface was taken to be 
203 K, to match the initial state in the sh ock experiments. Th e rocky core was modeled 



using an EOS fo r basalt, SESAME 7530 flBarnes fc Lvon 



19881 ). as was done previously 



( ISwift et al 



201lh . 




Fig. 2. — Principal isentrope of liquid ammonia (initial temperature 203 K) deduced from 
mechanical equation of state fitted to principal shock Hugoniot, showing sensitivity to as- 
sumed Griinseisen parameter. 
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The variations in F made a negligible difference to the mass-radius relations. At high 
masses, the mass-radius relation for pure ammonia asymptoted to a power-law behavior 
R = aMl^ with a = 1.4395 ± 0.0005 and (3 = 0.32889 ± 0.00004. For an incompressible 
material, 

M = ^vrrVo, (5) 

giving /3 = 1/3. The difference in the fitted value is small but significant; a is considerably 
less than the incompressible value of (|7rpo) • The mass-radius relation was also deduced 
using the SESAME EOS, which matched that from the Griineisen EOS up to 0.1 Mg, above 
which point the extrapolation beyond the bounds of the table gave unphysical behavior. 
(Figs [3] to El) 

The planetary radius for pure ammonia did not exhibit a maximum within the range 
of masses investigated. The range of compressions explored by the shock experiments was 
equivalent to the central pressure in bodies of pure ammonia up to around 2/3 Mg(1.5 Rjr;). 
However, the mass-radius relation is accurate for significantly larger bodies, because the 
mass and volume are dominated by matter at much lower pressures until the average density 
exceeds 1.5g/cm^ or so: approximately 4Mg, 2.5 Rg, and a core pressure of 100 GPa. The 
relation may be accurate for even larger bodies, but it has not been validated by EOS 
experiments. 

4. Conclusions 

The relation between shock and particle speeds in liquid ammonia appears linear to 
within the scatter in the data up to pressures of at least 39 GPa. A Griineisen mechanical 
equation of state was constructed using the principal Hugoniot of initial state zero pressure 
and 203 K as a reference, and estimating the Griineisen parameter F from the slope of the 
Hugoniot. Isentropes were calculated through the same state, the sensitivity to F rising 
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Fig. 3. — Mass-radius relation deduced from equation of state for liquid ammonia (surface 
temperature 203 K), also showing relation for incompressible material and least-squares fit 
to the relation, which is dominated by high pressure behavior. 
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- 11 - 




Fig. 5. — Mass-radius relations for differentiated ammonia-rock bodies. The percentages are 
the mass fraction of ammonia in the body. 
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with pressure. 

Mass-radius relations were calculated for self-gravitating bodies consisting of ammonia, 
and differentiated ammonia-rock mixtures. The mass-radius relations were insensitive 
to variations in F, indicating that the relations should be reliable for comparison to 
planetary measurements, to central pressures substantially above those reached in the shock 
experiments. 
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